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BACKGROUND OF THE INVENTION 
;g l . FIELD OF THE INVENTION 
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[02] The present invention relates to a text mining method and apparatus for 
ill extracting features of documents. In particular, the invention relates to a text mining 

y method and apparatus for extracting features of documents, wherein features are 

, = v. 

U extracted such that all mutually associated documents and terms are placed near each 

other in the feature space. Applications of the invention include document and/or web 
retrieval, associated term retrieval, document classification. 
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|l|§So 2. DESCRIPTION OF THE RELATED ART 
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[03] In text mining as a technology for squeezing desired knowledge or 
information by making analysis of text data, effective feature extraction of documents is 
an important task for efficiently performing document and/or web retrieval, associated 

- 1 - Docket No. 14321.31 



term retrieval, document classification and so on. As a typical document feature 
extracting method, a vector-space model as set out on page 3 1 3 of "Automatic Text 
Processing" (Addison- Wesley, 1989) is frequently used. 

[04] In the vector-space model, when terms selected as indices in the documents, 
namely index terms representing the contents of the documents, are t in number, a vector 
Vi is used respectively to correspond to an index term Tj to define a t-dimensional vector 
space. All vectors forming thus the defined vector space can be expressed as a linear 
!y combination of t in number of the vectors corresponding to t in number of the index 

1=^ terms. In this vector space, a document D r is expressed as follows: 

m = Z x (r V ; (1) 

5 / = 1 

\* j [05] In the foregoing expression (1), x lT active on Vj is the contribution of the index 

i 

S3 term Ti to the document D r and represents a feature of the document. The feature is an 

amount representing the term frequency of the index term in the document. A vector [x r i , 
>- x r 2, ... Xn]' of t x 1 (t rows and one column) becomes a feature vector of the document D r . 

jO — 

2 § 5 R 1 5 As the simplest case, when the index term Ti appears in the document D r , Xi r is set to 1 . 

g 8 < g * 5 . 

z § 1 oh § When the index term Tj does not appear in the document D r , Xj r is set to 0. In a more 

■7* CO O < Vi 2 

§ £ " § complicated case, as set forth in the foregoing publication on page 279 to 280, two 
quantities are used. These two quantities are a term frequency tf„ of the index term Tj in 
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the document D r and a document frequency dfj of documents containing the index term Ti 
in all documents registered in the document database. 

[06] For the group of documents consisting of d in number of documents, a t * d 
term-document matrix X can be defined as follows: 

X = [Xi, X2, X<j] 

[07] Here, a t-dimensional vector Xj = [xji, Xj2,..., xjj 1 expresses the feature vector of 
the document Dj, and '(dash) represents matrix inversion. 

[08] Fig. 1 is an illustration showing one example of documents, translated from 
Japanese sentences, registered in a document database, where "ronin" is a romanized 
word meaning students who, having failed a school entrance-exam of a particular 
academic year, are preparing for one next year. Fig. 2 is an illustration showing one 
example of a term-document matrix taking the Kanji (Chinese) characters appearing on 
the documents shown in Fig. 1 as index terms. Kanji terms- are underlined in Fig. 1. In 
Fig. 2, among a character string "let me know about" appearing in all of the documents 1 



2 > S § s to 3, the Kanji term "know" is checked off from the index terms. Fig. 3 is an illustration 

£ i ° 1 5 



p i & ° 1 h" showing one example of an actual input question, translated from Japanese, from a user, 

z 2 s o h w 

2 " c/5 O < v> ^ 
w H UJ < < 

1 o < 1 8 5 where Kanji terms are underlined. If the index terms of Fig. 2 are used to express the 



question, the question can be expressed with the term-document matrix shown in Fig. 4. 
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[09] 



In general, when the vector-space model is used, similarity sim (D r , D s ) of two 



documents D r and D s can be expressed as follows: 



^ X i> X is 



sim{D r ,D s ) = 



i = 1 



(2) 



Jl] 4 1 4 

V i = 1 / = 1 




[010] When the similarity of the question and each document of Fig. 1 is judged on 
the basis of the meaning of the question of Fig. 3, the question of Fig. 3 is the most similar 
to the document 3 of Fig. 1 . However, using the feature vectors as shown in Figs. 2 and 4, 
the similarity of each document of Fig. 1 and the question of Fig. 3 is respectively 
sim(document 1, question) = 0.5477, sim(document 2, question) = 0.5477, sim(document 
3, question) = 0.5477. In short, all have the same similarity. 

[011] As a solution for such a problem, a method called Latent Semantic Analysis 
(LSA) was proposed in "Journal of the American Society for Information Science" 1990, 
Vol. 41, No. 6, pp. 391 to 407. This method extracts latent meaning of the documents on 
the basis of co-occurrences of the terms and is significantly outstanding in terms of 
retrieving efficiency. Here, "co-occurrences of terms" represents a situation where the 
terms appear simultaneously in the same documents/statements. 

[012] The LSA extracts a latent semantic structure of the documents by performing 
singular value decomposition (SVD) for the term-document matrix. In the obtained 
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feature space, mutually associated documents and terms are located near each other. In a 
report placed in "Behavior Research Methods Instruments & Computers" (1991), Vol. 
23, No. 2, pp. 229 to 236, retrieval using the LSA indicates a result of 30% higher 
efficiency in comparison with the vector-space model. LSA will be explained hereinafter 
in more detail. 

[013] In LSA, at first, singular value decomposition is performed for the t x d 
term-document matrix X as set out below. 

X = T 0 SoDo' (3) 
[014] Here, To represents an orthogonal matrix of t * m, So represents a square 
diagonal matrix of m x m with taking m in number of the singular values as the diagonal 
elements and setting 0 to the other elements. Do' represents an orthogonal matrix of m x 
d. In addition, let us assume that 0 > d > t, and arrange the orthogonal elements of So in 
descending order. 

[015] Furthermore, in LSA, with respect to the feature vector x q of t x 1 of a 



1 > § 1 s document D q , the following conversion is performed to derive a LSA feature vector y q of 



gdSjSSfc: n xl: 



Z O ct o f- U3 
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I|<1sh y^s-^x, ( 4) 
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[016] Here, S is a square diagonal matrix of n x n taking the first to (n)th of the 
diagonal elements of So, and T is a matrix of t x n drawing the first to (n)th columns of To. 
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[017] As an example, results of singular value decomposition of the term-document 
matrix shown in Fig. 2 are given below. The matrices To, So and Do are expressed as 



follows: 



0.1787 
0.1787 
0.1787 
0.4314 
0.4314 
0.1787 
0.1787 
0.4314 
0.4314 
0.1787 
0.2527 



S n = 



-0.3162 
-0.3162 
-0.3162 
-0.3162 
-0.3162 
0.3162 
0.3162 
0.3162 
0.3162 
0.3162 
0.0000 



2.7979 0 
0 2.2361 
0 0 



0.3393 
0.3393 
0.3393 
-0.1405 
-0.1405 
0.3393 
0.3393 
-0.1405 
-0.1405 
0.3393 
0.4798 

0 
0 

1.4736 



Do = 



0.5000 -0.7071 0.5000 
0.5000 0.7071 0.5000 
0.7071 0.0000 -0.7071 
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[018] Let us assume that the dimension t of the LSA feature vectors is 2 and 
applying the foregoing expression (4) to each feature vector of the term-document matrix 
in Fig. 2. Then, the LSA feature vectors of the documents 1, 2 and 3 are respectively 
[0.5000, -0.7071]', [0.5000, 0.7071]' and [0.7071, 0.0000]'. In addition, applying the 
foregoing expression (4) to the feature vector of Fig. 4, the LSA feature vector of the 



question from the user becomes [0.6542, 0]'. 
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[019] Applying the foregoing expression (2) to the LSA feature vectors obtained as 
set forth above, the similarity of the question of Fig. 3 and each document of Fig. 1 , 
become respectively, sim(document 1, question) = 0.5774, sim(document 2, question) = 
0.5774, and sim(document 3, question) = 1.0000. Thus, a result that the document 3 has 
the highest similarity to the question can be obtained. Considering a help system 
application or the like utilizing computer networks, an answer statement of the document 
3 registered in the document database will be returned to the user who asked the question 
of Fig. 3. 

[020] For singular value decomposition, an algorithm proposed in "Matrix 
Computations", The Johns Hopkins University Press, 1996, pp. 455 to 457, is frequently 
used. In the report of "Journal of the American Society for Information Science" set forth 
above, there is a statement that the value of the number of rows (or columns) n of the 
square matrix S is preferably about 50 to 150. In addition, in the foregoing report of 
"Behavior Research Methods, Instruments, & Computers", it has been indicated that 
better efficiency can be attained by pre-processing using the term frequency or document 
frequency instead of simply setting each element of the feature vector to 0 or 1 before 
performing LSA. 

[021] However, in the algorithm for singular value decomposition proposed in the 
foregoing "Matrix Computations", memory space in the order of the square of the number 
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of index terms t (t 2 ) is required at the minimum. This is because a matrix of t * t is 
utilized for bidiagonalization of a matrix in the process of calculation of basis vectors 
spanning a feature space from a given term-document matrix. The prior art is therefore 
not applicable to document database holding a very large number of terms and data. 
Furthermore, the prior art requires complicated operations of matrices irrespective of the 
number of data. 



in 
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SUMMARY OF THE INVENTION 
[022] The present invention has been worked out in view of the problems set forth 
above. In a first aspect of the present invention, a text mining method is provided for 
extracting features of documents using a term-document matrix consisting of vectors 
corresponding to index terms representing the contents of the documents. In the 
term-document matrix, contributions of the index terms to each document act on 
respective elements of the term-document matrix. The method comprises: 
1 J a basis vector calculating step of calculating a basis vector spanning a 

1^ feature space, in which mutually associated documents and terms are located in 

jL§ proximity with each other, based on a steepest descent method minimizing a cost; 

a feature extracting step of calculating a parameter for normalizing 
^ features using the term-document matrix and the basis vector, and extracting the 

features on the basis of the parameter; and 

a term-document matrix updating step of updating the term-document 



1 5 1 3 3 matrix to a difference between the term-document matrix, to which the basis 

ot 5 -J H £3 < 
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vector is not applied, and the term-document matrix, to which the basis vector is 



3 < I s h applied. 



[023] In a second aspect of the present invention, a text mining method is provided 
for extracting features of documents wherein the cost is defined as a second-order cost of 
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the difference between the term-document matrix, to which the basis vector is not 
applied, and the term-document matrix, to which the basis vector is applied. 
[024] In a third aspect of the present invention, a text mining method is provided for 
extracting features of documents wherein the basis vector calculating step comprises: 

an initializing step of initializing a value of the basis vector; 

a basis vector updating step of updating the value of the basis vector; 

a variation degree calculating step of calculating a variation degree of the 
value of the basis vector; 

a judging step of making a judgment whether a repetition process is to be 
terminated or not using the variation of the basis vector; and 

a counting step of counting the number of times of the repetition process. 
[025] In a fourth aspect of the present invention, a text mining method is provided 
for extracting features of documents wherein the basis vector updating step updates the 
basis vector using a current value of the basis vector, the term-document matrix and an 
updating ratio controlling the updating degree of the basis vector. 
[026] In a fifth aspect of the present invention, a text mining method is provided for 
extracting features of documents wherein when all basis vectors and normalizing 
parameters required in extracting the features have been already obtained, the calculation 
of the normalizing parameters in the basis vector calculating step and the execution of the 
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feature extracting step are omitted. In addition, the feature extracting step extracts the 
features using the basis vectors and the normalizing parameters that have been already 
obtained. 

[027] In a sixth aspect of the present invention, a text mining apparatus is provided 
for extracting features of documents using a term-document matrix consisting of vectors 
corresponding to index terms representing the contents of the documents. In the 
term-document matrix, contributions of the index terms to each document act on 
;?? respective elements of the term-document matrix. The apparatus comprises: 

!=* basis vector calculating means for calculating a basis vector spanning a 

\u feature space, in which mutually associated documents and terms are located in 

proximity with each other, based on a steepest descent method minimizing a cost; 
t ^ feature extracting means for calculating a parameter for normalizing 

features using the term-document matrix and the basis vector, and extracting the 
features on the basis of the parameter; and 

S§ = 

" 3 £ I s! s term-document matrix updating means for updating the term-document 

* 9 i o 2 3 

Otf 5 J H w < 
o g < p x 

§ ^ > o d £ matrix to a difference between the term-document matrix, to which the basis 

r 5 S u o n 

I S < I g h vector is not applied, and the term-document matrix, to which the basis vector is 



applied. 
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[028] In a seventh aspect of the present invention, a text mining apparatus is 

provided for extracting features wherein the cost is defined as a second-order cost of the 

difference between the term-document matrix, to which the basis vector is not applied, 

and the term-document matrix, to which the basis vector is applied. 

[029] In an eighth aspect of the present invention, a text mining apparatus is 

provided for extracting features of documents wherein the basis vector calculating means 

comprises: 

^ initializing means for initializing a value of the basis vector; 

\^ basis vector updating means for updating the value of the basis vector; 

|y variation degree calculating means for calculating a variation degree of 

\% the value of the basis vector; 

^ judging means for making a judgment whether a repetition process is to 

be. terminated or not using the variation of the basis vector; and 

counting means for counting the number of times of the repetition 

> - 

So = 

£ < £ g 3 « process. 

ot h J t- S <. 

g ^ > o 5 £ [030] In a ninth aspect of the present invention, a text mining apparatus is provided 

g 5 5 3 8 u 

2 co P < co *d 
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1 1 < i s h for extracting features of documents wherein the basis vector updating means updates the 
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basis vector using a current value of the basis vector, the term-document matrix and an 
updating ratio controlling the updating degree of the basis vector. 
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[031] 



In a tenth aspect of the present invention, a text mining apparatus is provided 




for extracting features of documents wherein when all basis vectors and normalizing 
parameters required in extracting the features have been already obtained, the calculation 
of the normalizing parameters in the basis vector calculating means and the execution of 
the feature extracting means are omitted. In addition, the feature extracting means 
extracts the features using the basis vectors and the normalizing parameters that have 
been already obtained. 

[032] In an eleventh aspect of the present invention, a computer program product is 
provided for being executed in a text mining apparatus for extracting features of 
documents using a term-document matrix consisting of vectors corresponding to index 
terms representing the contents of the documents. In the term-document matrix, 
contributions of the index terms act on respective elements of the term-document matrix. 
The computer program product comprises: ~ 

a basis vector calculating step of calculating a basis vector spanning a 

feature space, in which mutually associated documents and terms are located in 

proximity with each other, based on a steepest descent method minimizing a cost; 
a feature extracting step of calculating a parameter for normalizing 

features using the term-document matrix and the basis vector, and extracting the 

features on the basis of the parameter; and 
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a term-document matrix updating step of updating the term-document 
matrix to a difference between the term-document matrix, to which the basis 
vector is not applied, and the term-document matrix, to which the basis vector is 
applied. 

[033] The feature extracting apparatus disclosed in this specification is constructed 
by defining a cost as a second-order function of a difference between the term-document 
matrix, to which the basis vector is not applied, and the term-document matrix, to which 
the basis vector is applied. The apparatus merely requires the following means: 

a. basis vector calculating means for calculating a basis vector 
applying a steepest descent method to the cost; 

b. feature extracting means for calculating a parameter for 
normalizing features using the term-document matrix and the basis vector, and 
extracting the features on the basis of the parameter; 

c. term-document matrix updating means for updating the 



P 3 1 ° ^ which the basis vector is not applied, and the term-document matrix, to which the 



basis vector is applied, in order to prevent redundant extraction of features; and 

d. feature extraction control means for controlling execution of 
respective means. 
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[034] The basis vector calculating means repeats calculation on the basis of the input 
term-document matrix to finally derive one basis vector. The repetition process is 
terminated when the variation degree of the basis vector becomes less than or equal to a 
predetermined reference value. 

[035] The feature extracting means calculates a parameter for normalizing the 
features on the basis of the input basis vector and the term-document matrix, and extracts 
a feature for each document. 

[036] The term-document matrix updating means updates the term-document matrix 
on the basis of the input basis vector. 

[037] The feature extraction control means repeats execution of each means until the 
number of the features defined by the user is satisfied. When the basis vectors and the 
normalizing parameters have been already calculated, execution of the basis vector 
calculating means and calculation of the normalizing parameters in the feature extracting 
means are omitted. Then, the feature extraction can be performed with the construction 
incorporating the already obtained basis vectors and the normalizing parameters. 
[038] According to the present invention, a text mining method for extracting 
features of documents using a term-document matrix consisting of vectors corresponding 
to index terms representing the contents of the documents, wherein contributions of the 
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index terms act on respective elements of the term-document matrix, comprises the 
following steps: 

i. a basis vector calculating step of calculating a basis vector 
spanning a feature space, wherein mutually associated documents and terms are 
located in proximity with each other based on a steepest descent method 
minimizing a cost; 

ii. a feature extracting step of calculating a parameter for normalizing 
features using the term-document matrix and the basis vector, and extracting the 
features on the basis of the parameter; 

iii. a term-document matrix updating step of updating the 
term-document matrix to a difference between the term-document matrix, to 
which the basis vector is not applied, and the term-document matrix, to which the 
basis vector is applied; and 

iv. a feature extraction control step of controlling execution of 
respective steps. 

[039] Therefore, concerning feature extraction of documents in text mining, the 
features having the same nature as those obtained by LSA can be extracted with smaller 
memory space than the apparatus or method executing LSA. In addition, specific 
software or hardware for extracting the features can be easily implemented. 
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[040] The above and other objects, characteristics and advantages of different 
embodiments of the present invention will become more apparent from the following 
descriptions of embodiments thereof taken in conjunction with the accompanying 
drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[041] Fig. 1 is an illustration showing one example of documents registered in a 
document database; 

[042] Fig. 2 is an illustration showing one example of a term-document matrix with 

taking Kanji terms appearing in the documents shown in Fig. 1 as index terms; 

[043] Fig. 3 is an illustration showing one example of a question actually input by a 

user; 

[044] Fig. 4 is an illustration showing a term-document matrix obtained from the 
question in Fig. 3; 

[045] Fig. 5 is an illustration showing one embodiment of a feature extracting 
apparatus according to the present invention; 

[046] Fig. 6 is an illustration showing one example of a hardware construction for 
implementing the present invention; 

[047] Fig. 7 is an illustration showing a structure of a term-document matrix data 



J*I1 S file; 



g 3 > o 5 £ [048] Fig. 8 is an illustration showing a structure of a basis vector data file, wherein 



v/ o < 1 s i- the calculated basis vectors are stored: 



[049] Fig. 9 is an illustration showing a structure of a feature data file; 
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[050] Fig. 10 is an illustration showing a structure of a normalizing parameter data 
file; 

[051] Fig. 1 1 is a flowchart showing calculation of a basis vector in basis vector 
calculating means; and 

[052] Fig. 12 is an illustration showing one example of an automatic document 
classifying system employing one embodiment of the feature extracting apparatus 
according to the present invention. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[053] Fig. 5 is an illustration showing one embodiment of a feature extracting 
apparatus according to the present invention. As shown in Fig. 5, feature extraction 
control means 200 has term-document matrix updating means 210, basis vector 
calculating means 220, feature extracting means 230. The reference numeral 100 denotes 
a term-document matrix data file, 300 denotes a basis vector data file, 400 denotes a 
feature data file, 450 denotes a normalizing parameter data file. In the term-document 
matrix data file 100, a term-document matrix of collected document data is stored. The 
term-document matrix updating means 210 reads the term-document matrix from the 
term-document matrix data file 1 00, and transfers the read term-document matrix to the 
basis vector calculating means 220 and the feature extracting means 230 without updating 
the term-document matrix, in a first iteration process. 

[054] In a second and subsequent iteration processes, the terms-document matrix is 
updated on the basis of the basis vector transferred from the basis vector calculating 
means 220. The result of updating is transferred to the basis vector calculating means 220 
and the feature extracting means 230. The basis vector calculating means 220 calculates 
one basis vector through a repetition process based on the term-document matrix 
transferred from the term-document matrix updating means 210. Then, the degree of 
variation of the basis vector in respective repetition is monitored for terminating the 
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repetition process when the degree of variation becomes less than or equal to a 
predetermined reference value. The basis vector calculating means 220 stores the 
calculated basis vector in the basis vector data file 300 and in conjunction therewith, 
transfers the calculated basis vector to the term-document matrix updating means 210 and 
the feature extraction means 230. The feature extracting means 230 extracts one feature 
of each document on the basis of the term-document matrix transferred from the 
term-document matrix updating means 210 and the basis vector transferred from the basis 
:J vector calculating means 220. The result is stored in the feature data file 400, and also the 

\^ parameter for normalizing the features is stored in the normalizing parameter data file 

[055] Execution of the term-document matrix updating means 210, the basis vector 
|^ calculating means 220 and the feature extracting means 230 is taken as one iteration 

process. Number of times of iteration processes will be indicated by suffix i, and number 
of features designated by the user is indicated by suffix n. The feature extraction control 
s 2 £ 1 3 s means 200 repeats the process until a condition, i = n is satisfied. On the other hand, in a 
p ^ 1 2 o t case w ^ ere a ^ °f required basis vectors and the required normalizing parameters have 



z 



r* w3 O < co U 



<ujh 0 uj 



§< 



^^Is^ already been obtained, execution of the basis vector calculation means 220 and 
calculation of the normalizing parameters in the feature extracting means 230 may be 
omitted. Therefore, in such a case, the feature extraction control means 200 may be 
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constructed with the term-document matrix updating means 210 incorporating the known 
basis vectors and normalizing parameters, and with the feature extracting means 230. 
[056] Fig. 6 is an illustration showing one example of a hardware construction for 
implementing the present invention. As shown in Fig. 6, the feature extracting apparatus 
includes the following components: 

a central processing unit (CPU) 10 performing control for the overall 
apparatus, 

a memory 20 for storing the program and providing a temporary data 
storage region required for executing the program, 

a keyboard 30 for inputting data, and 

a display 40 for generating a display screen. 
[057] The programs to be executed by the feature extraction control means 200, the 
term-document matrix data file 100, the basis vector data file 300, the feature data file 
400, and the normalizing parameter data file 450 are stored in the memory 20. 
[058] By taking this construction, the feature extraction is performed by CPU 10 
receiving the command from the user through the keyboard 30, a mouse pointing a 
desired position on the display 40, or the like. It should be noted that, in the example 
shown in Fig. 5, the feature extraction control means 200 has a stand-alone construction. 
However, the feature extraction control means 200 may be built-in other systems. 
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[059] Fig. 7 is an illustration showing a structure of the term-document matrix data 
file. In Fig. 7, the reference numerals 101-1, 101-2, 101-d correspond to t-dimensional 
term-document data A consisting of d in number of data. Here, X = [xj, X2, ... xj, Xj = 
[xji, Xj 2 , ... Xj t ]' are defined to express the term-document data A with a t * d matrix X. 
[060] Fig. 8 is an illustration showing a structure of the basis vector data file storing 
the calculated basis vectors. In Fig. 8, the reference numerals 301-1, 301-2, ... 301-n 
correspond to t-dimensional basis vector data B consisting of n in number of data. The 
(i)th element 301-i corresponds to an output value of the basis vector calculating means 
220 in the (i)th iteration process in Fig. 5. In the following disclosure, this element is 
expressed by a t x 1 column vector Wj = [wn, Wj2, Wj t ]\ 

[061] Fig. 9 is an illustration showing a structure of the feature data file. In Fig. 9, 
the reference numerals 401-1, 401-2, 401-n correspond to d-dimensional feature data 
C consisting of n in number of data. - The (i)th element 401— i corresponds to an output 
value of the feature by the feature extraction means 230 in the (i)th iteration process in 



s 2 > I § s Fig- 5. This element is expressed by an 1 x d row vector y\ = [yn, y-a, yid]. 



p i. & ° 3 h [062] Fig. 1 0 is an illustration showing a structure of the normalizing parameter data 
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| fe < 1 1 h file. In Fig. 1 0, the reference numerals 45 1 - 1 , 45 1 -2, 45 1 -n correspond to normalizing 
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parameter data D consisting of n in number of data. The (i)th element 45 1 -i corresponds 
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to an output value of the normalizing parameter by the feature extracting means 230 in the 
(i)th iteration process in Fig. 5. 

[063] Using the foregoing definitions, an implementation of feature extraction in the 
shown embodiment will be explained. The term-document matrix updating means 210 
reads out X from the term-document matrix data file 100 only when i = 1, namely in the 
first iteration process, to store in a t * d matrix E without performing any arithmetic 
operation. Accordingly, E = [e u e 2 , e d ], ej = [eji, e j2 , e jt ]' = [x jl5 x j2 , x jt ]\ In order 
to prevent redundant extraction of the features extracted in the preceding iteration 
processes, E is updated in the (i)th iteration using the current value and the basis vector 
calculated in the immediately preceding iteration process. The result of updating is 
transferred to the basis vector calculating means 220. A value of E in the (i)th iteration, 
E(i), will be expressed by the following expression (5): 



[064] Here, E(i) = [ej(i), e 2 (i), e d (i)], each element e/i) of E(i) is defined by ej(i) = 
[eji(i), ej 2 (i), ej t (i)]'. Namely, when i > 2, the term-document matrix is updated to a 
difference derived by subtracting the term-document matrix, to which the basis vector is 
applied, from the term-document matrix, to which the basis vector is not applied. 
[065] Fig. 1 1 is a flowchart showing calculation of the basis vector in the basis 
vector calculating means. In Fig. 1 1, a value of Wj in the (k)th repetition is expressed by 




X, for / = 1 

E(/-l) - Wj.^w'j.! E(/-l)), otherwise 



(5) 
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Wi(k) = [wii(k), w^Ck), w it (k)]'. At first, at step S500, the suffix k is initialized to 1. 
Subsequently, the process is advanced to step S510 to initialize respective element of 
Wj(l) with an arbitrary value between -C to C. Here, the value of C may be a positive 
small value, such as C = 0.01. At step S520, in order to calculate the basis vector 
spanning a feature space where mutually associated documents and terms are located in 
proximity with each other, a second-order cost expressed by the following expression (6) 
is provided. 

^1 £(e,„(0- w /( .y„) 2 (6) 

m = l / = ) 

[066] Here, "terms are placed in proximity" means that the positions of the terms are 
close with each other within a feature space, and "documents are placed in proximity" 
means that the positions of terms included in respective documents are close in the 
feature space. On the other hand, a cost means an object to be minimized. In the shown 
embodiment, the cost is defined as a second-order function of the difference between the 
term-document matrix, to which the basis vector is not applied, and the term-document 
matrix, to which the basis vector is applied, as expressed by the expression (6). Here, 



is the (m)th element of a 1 x d vector y i which is defined as follows: 

y, =[yi,,yi 2 > 5 y id ] = w/eco (7) 
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[067] For the cost, the steepest descent method is applied to update the value of w\ as 
expressed by the following expression (8). 

w.tf + l) = w,(*) + ^(E(i) - z,{kY (8) 

a 

[068] Here, ix i(k) is an update ratio controlling the degree of updating in the (k)th 
repetition, which is initialized by a positive small value when k is 1 , such as \x \{X) = 0.1 . 
Every time of increment of k, the value is decreased gradually. In the alternative, it is also 
possible to set the value at a constant value irrespective of k. On the other hand, Zj(k) is 
defined as follows: 

Zi (k) = Wi(k)'E(i) (9) 
[069] At step S530, 5 i(k) indicating the degree of variation of Wj is derived as 
follows: 

= " W ^ W)I (10) 

[070] At step S540, a judgment is made whether the process is to be terminated or 
not on the basis of the value of 5 i(k). As a result of the judgment, if termination is 
determined, the process is advanced to step S560, and otherwise, the process is advanced 
to step S550. Here, in Fig. 11, j3 j is a positive small value, such as j3 \ = 1 * 10~ 6 . 
[071] At step S550, the value of the counter k is incremented by 1 . Then, the process 
is returned to step S520. At step S560, w\ is stored as the (i)th data of the basis vector data 
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file 300. At the same time, Wj is transferred to the term-document matrix updating means 
2 1 0 and the feature extracting means 230. In the feature extracting means 230, the feature 
and the normalizing parameter p\ are calculated in the following manner. 

y, =y,//>, (ii) 

[072] Here, p x is defined as follows: 



d 

2 



Pi = JZ Vu 

V Jml (12) 
[073] The feature y\ and the normalizing parameter p\ are stored respectively in the 

feature data file 400 and the normalizing parameter data file 450 as the (i)th data. 

[074] Fig. 12 is an illustration showing one example of an automatic document 

classifying system employing the shown embodiment of the feature extracting apparatus. 

In Fig. 12, the reference numeral 601 denotes term-document matrix calculating means, 

602 denotes classifying means. The classifying means 602 may be implemented by a 

method disclosed in "Journal of Intelligent and Fuzzy Systems", published on 1993, Vol. 

l,No. 1, Pages 1 to 25. 



1 8 5 g § 5 [075] The document data stored in the document database E is taken in the automatic 

4J s f- u P b 

z § 1 1 h w document classifying system 600. In the automatic document classifying system 600, a 

Zr Y\ t W <i < 

|£ < term-document matrix is derived in the term-document matrix calculating means 601. 



The result of calculation of the term-document matrix is transferred to the feature 
extraction control means 200. The feature extraction control means 200 extracts the 



features from the received term-document matrix. The extracted result is output to the 
classifying means 602. In the classifying means 602, the result of classification is output 
on the basis of the input features. 

[076] To evaluate the present invention, feature extraction of actual document data 
related to an entrance examination system was performed. It has been confirmed that the 
present invention could extract the features of the same nature as those extracted using 
the conventional LSA. 

[077] Next, concerning the size of the memory space, in a typically practical case 
where the number of terms t is significantly greater than number of documents d (t » d), 
the conventional LSA requires in the order of t of the memory size, the present invention 
merely requires the memory size in the order of t-d for calculating respective basis 
vectors. Furthermore, in order to realize the prior art, a complicated matrix operation 
apparatus is required. _ The system according to the invention, however, can be easily 
realized with an apparatus that performs simple arithmetic operations. Namely, 
according to the present invention, the LSA feature extraction can be performed using a 
smaller memory space and a simpler program. In addition, this simple program may be 
loaded in a digital signal processor (DSP). Therefore, a specific chip for feature 
extraction can be produced easily. 
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[078] Hereinafter, the results of respective means executing the shown embodiment 
of the feature extracting apparatus for the documents of Fig. 1 and the question of Fig. 3 



will be shown. 



[079] A. Documents of Fig. 1 



[080] First, let X denote the term-document matrix of Fig. 2. 



[081] I. First Iteration in Feature Extraction Control Means 200 (i = 1) 

[082] According to the foregoing expression (5), the term-document matrix 

updating means 210 outputs E(l) expressed by the following expression to the basis 



vector calculating means 220 and the feature extracting means 230. 



E(l) = 



1 0 0 

1 0 0 

1 0 0 

1 0 1 

1 0 1 

0 1 0 

0 1.0 

0 1 1 

0 1 1 

0 1 0 

0 0 1 



[083] In the basis vector calculating means 220, initialization is performed with 
setting the basis vector wi(l) at [0.0100, -0.0100, 0.0100, -0.0100, 0.0100, -0.0100, 
0.0100, -0.0100, 0.0100, -0.0100, 0.0 100]', n i at a fixed value 0.1, 0 i at 1 x 10' 6 . The 
calculation shown in Fig. 1 1 is performed for a hundred thirty-two times. Then, the basis 
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vector w, = [0.1787, 0.1787, 0.1787, 0.4314, 0.4314, 0.1787, 0.1787, 0.4314, 0.4314, 
0.1787, 0.2527]'is stored in the basis vector data file 300, and output to the feature 
extracting means 230 and the term-document matrix updating means 210. 
[084] * First Repetition in Basis Vector Calculating Means 220 (k = 1) 
[085] From the foregoing expression (8), 

w,(2) = [0.0103 -0.0097, 0.0103, 0.0093, 0.0107, -0.0103, 0.0097, -0.0100, 0.0100, 
-0.0103,0.0103]' 

wi(2) - wi(l) = 10' 3 x [0.3332, 0.3334, 0.3332, 0.6668, 0.6666, -0.3332, -0.3334, 0.000-1, 
-0.0001,-0.3332,0.3332]' 

5 i(l) = 0.0103 

[086] * Second Repetition in Basis Vector Calculating Means 220 (k = 2) 
[087] From the foregoing expression (8), 

w,(3) = [0.0107, -0.0093, 0.0107, -0.0085, 0.0115, -0.0107, 0.0093, -0.0100, 0.0100, 
-0.0107, 0.0107]' 

w,(3) - w,(2) = 10' 3 x [0.41 10, 0.41 12, 0.41 10, 0.8001, 0.7998, -0.3665, -0.3668, 0.0224, 
0.0221,-0.3665,0.3887]' 

5 ,(2) = 0.0015 
[088] — syncopated — 
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[089] * A Hundred Thirty-Second Repetition in Basis Vector Calculating Means 
220 (k= 132) 

[090] From the foregoing expression (8), 

w,(133) = [0.1787, 0.1787, 0.1787, 0.4314, 0.4314, 0.1787, 0.1787, 0.4314, 0.4314, 
0.1787, 0.2527]' 

w,(133) - wi(132) = 10' 6 x [-0.3020, -0.3020, -0.3020, -0.3020, -0.3020, 0.3020, 0.3020, 
0.3020, 0.3020, 0.3020, 0.0000]' 



5 ,(132) = 9.5500 x 10 



r7 



[091] 



In the feature extracting means 230, the operations shown in the expressions 



(11) and (12) are performed for outputting: 



y, =[0.5000, 0.5000, 0.7071] 



and 



pi =2.7979 



to the feature data file 400 and the normalizing parameter data file 450. 




expression (5), E(2) expressed as follows is output to the basis vector calculating means 



[092] 



[093] 



In the term-document matrix updating means 210, from the foregoing 



II. Second Iteration in Feature Extraction Control Means 200 (i = 2) 



220 and the feature extracting means 230: 
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E(2) = 



0 7^00 


-fi 7500 


-0 3536 


0 7S00 


-0 9S00 


-0 1S16 

U.J JJU 


0 7^00 


— 0 9S00 


— 0 1S16 

v. J JJU 


0 1Q64 


— 0 6016 


0 1464 




— v.DUjU 


0 1464 
W. 1 tot 


— o 9snn 


0 7S00 


_0 1<\1S 

— U.JJJJ 


-0.2500 


0.7500 


-0.3535 


-0.6036 


0.3965 


0.1465 


-0.6036 


0.3965 


0.1465 


-0.2500 


0.7500 


-0.3535 


-0.3536 


-0.3535 


0.5000 



[094] In the basis vector calculating means 220, initialization is performed with 
setting the basis vector w 2 (l) at [0.0100, -0.0100, 0.0100, -0.0100, 0.0100, -0.0100, 
0.0100,-0.0100,0.0100,-0.0100,0.0100]', // 2 at a fixed value 0.1, 0 2 atl x 10" 6 . The 
calculation shown in Fig. 1 1 is performed for a hundred nineteen times. Then, the basis 
vector w 2 = [0.3162, 0.3162, 0.3162, 0.3162, 0.3162, -0.3162, -0.3162, -0.3162, -0.3162, 
-0.3162, 0.0000]' is stored in the basis vector data file 300, and output to the feature 
extracting means 230 and the term-document matrix updating means 210. 
[095] * First Repetition in Basis Vector Calculating Means 220 (k = 1) 
[096] From the foregoing expression (8), 

w 2 (2) = [0.0102, -0.0098, 0.0102, -0.0096, 0.0104, -0.0105, 0.0095, -0.0103, 0.0097, 
-0.0105,0.0102]' 

w 2 (2) - w 2 (l) = 10" 3 x [0.2154, 0.2156, 0.2154, 0.3822, 0.3821, -0.451 1, -0.4513, -0.2844, 
-0.2846, -0.4511,0.1666]' 
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5 2 (1) = 0.001 1 



[097] * Second Repetition in Basis Vvector Calculating Means 220 (k = 2) 



[098] From the foregoing expression (8), 



w 2 (3) = [0.0105, -0.0095, 0.0105, -0.0092, 0.0108, -0.0110, 0.0090, -0.0106, 0.0094, 



-0.0110, 0.0103]' 



w 2 (3) - w 2 (2)= 10" 3 x [0.2624,0.2626, 0.2624, 0.4413, 0.441 1, -0.5152, -0.5154, -0.3364, 



So 

<6 O < 
otf 5 J 

88< 

^ co H 
< UJ f- 

^ < 

si 

5< 



UJ j 00 

H UJ < 

uj o n 

j w O 
O H UJ 

w < 5 

CO 



-0.3366, -0.5152, 0.1786]' 



5 2 (2) = 0.0013 



[099] ~ syncopated — 



[0100] * A Hundred Nineteenth Repetition in Basis Vector Calculating Means 220 (k 



119) 



[0101] From the foregoing expression (8), 



w 2 (120) = [0.3162, 0.3162, 0.3162, 0.3162, 0.3162, -0.3162, -0:3162, -0.3162, -0.3162, 



0.0000]' 



w 2 (120) - w 2 (119) = 10- 6 x [0.3327, 0.3333, 0.3327, -0.1375, -0.1381, 0.3332, 0.3326, 



-0.1377, -0.1383, 0.3332, -0.4712]' 



-7 



8 2 (1 19) = 9.8141 x 10 



[0102] In the f ea i ure extracting means 230, the operations shown in the expressions (1 1) 



and (12) are performed for outputting: 
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y 2 =[0.7071, -0.7071, -0.0000] 

and 

p 2 = 2.2361 

to the feature data file 400 and the normalizing parameter data file 450. 
[0103] From the results set forth above, the feature vectors of the documents 1 , 2 and 
3 in Fig. 1 are respectively [0.5000, 0.7071]', [0.5000, -0.7071], [0.7071, -0.0000]. 
Comparing these with the features of the LSA of respective documents shown in the 
explanation of the prior art, the second element of each vector is of opposite sign but has 
the same absolute value. Accordingly, concerning calculation of similarity in the 
expression (2), they have the same nature as the features of LSA. 
[0104] B. Question of Fig. 3 

[0105] Here, let us use the basis vectors stored in the basis vector data file 300 and the 
normalizing parameters stored in the normalizing parameter data file 450 during 
extraction of the features of the documents of Fig. 1. Thereby, execution of the basis 
vector calculating means 220 and calculation of the normalizing parameter in the feature 
extracting means 230 are omitted. Let X denote the term-document matrix of Fig. 4. 
[0106] I. First Iteration in Feature Extracting Means 200 (i = 1) 
[0107] In the term-document matrix updating means 210, E(l) expressed as follows 
from the foregoing expression (5) is output to the feature extracting means 230. 
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# 



0 



0 



E(l) = 



0 



0 



0 



[0108] In the feature extracting means 230, arithmetic operation according to the 
foregoing expressions (11) and (12) is performed using the feature vector wi and the 
normalizing parameter p\ obtained upon extraction of the features of the documents of 
Fig. 1 to output 

y, = [0.6542] 

to the feature data file 400. 
[0109]. II. Second Iteration in Feature Extraction Control Means 200 (i = 2) 
[0110] In the term-document matrix updating means 210, using the feature vector wj 
obtained upon performing feature extraction of the documents shown in Fig. 1 , from the 
foregoing equation (5), E(2) expressed as follows is output to the feature extracting 
means 230. 
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• 



E(2) = 



-0.3271 
-0.3271 
0.6729 
0.2103 
0.2103 
0.6729 
0.6729 
-0.7897 
0.2103 
-0.3271 
-0.4626 



[0111] In the feature extracting means 230, arithmetic operation according to the 
foregoing expressions (11) and (12) is performed using the feature vector W2 and the 
normalizing parameter pi obtained upon extraction of the features of the documents of 



Fig. 1 to output 



y 2 = [-0.0000] 



to the feature data file 400. 



[0112] From the result set forth above, the feature vector of the question of Fig. 3 
becomes [0.6542, -0.0000]', comparing the value explained in the prior art, the second 



element has the same absolute value. 



[0113] The present invention has been described in detail with respect to preferred 
embodiments. It will now be apparent from the foregoing to those skilled in the art that 



changes and modifications may be made without departing from the invention in its 
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broader aspect. It is the intention, therefore, in the apparent claims to cover all such 
changes and modifications as fall within the true spirit of the invention. 
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